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This  paper  describes  the  power  generation  behavior  of  a  thermoelectric  generator  (TEG)  exposed  to 
a  transient  heat  source  on  the  hot  side  and  natural  convection  on  the  cold  side.  The  simulation  situation 
is  typical  in  energy  harvesting  applications.  Modeling  thermoelectric  generators  (TEGs)  under  these 
conditions  is  complicated  compared  to  thermoelectric  coolers  because  of  the  non-linearities  and  the 
unknown  electric  currents  in  a  closed-loop  circuit.  A  transient  thermoelectric  model  that  includes 
Seebeck,  Peltier,  Thomson,  and  Joule  effects  is  solved  using  finite-difference  techniques  and  the  power 
generated  from  a  TEG  is  simulated.  Using  open-circuit  experiments  were  used  to  establish  key  param¬ 
eters  governing  the  thermal  behavior  and  thermoelectric  coupling.  Experiments  with  closed-circuits  and 
load  resistors  were  used  to  validate  the  model.  The  results  show  that  inclusion  of  Thomson  effect  plays 
a  significant  role  in  accurately  predicting  the  power  generated  by  the  device. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  generators  (TEGs),  or  Seebeck  elements,  enable 
conversion  of  thermal  energy  into  electricity  without  any  moving 
parts.  The  Seebeck  elements  can  be  used  for  temperature  sensing  s 
and  power  generation  applications  [1].  TEGs  generate  electric 
current  in  a  closed-circuit  across  a  load  when  a  temperature 
gradient  develops  between  the  two  ends  of  the  device  [2-4].  The 
coupling  effects  governing  the  conversion  are  the  Seebeck,  Peltier, 
Thomson,  and  Joule  effects. 

Understanding  the  transient  behavior  of  TEG  is  important  for 
optimizing  energy  harvesting  from  waste  heat,  when  one  junction  is 
exposed  to  an  unsteady  heat  source  and  the  other  junction  is  sub¬ 
jected  to  natural  convection  at  ambient  temperature.  Several  models 
have  been  developed  for  such  configurations  [5—7],  but  most  models 
assume  either  an  unlimited  heat  source  or  steady  state  operational 
scenarios.  Practical  applications  such  as  parasitic  energy  harvesting 
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from  waste  heat  sources  generally  have  unsteady  heat  flux  and/or 
temperature  conditions  on  the  hot  side  of  the  TEG  and  pose  a  fully 
coupled  thermoelectric  problem  [8].  Simulation  of  TEGs  operating  in 
transient  mode  is  more  challenging  than  the  thermoelectric  coolers 
(TECs)  because  both  the  temperature  field  and  the  electric  currents 
vary  with  time.  The  problem  is  more  simplified  for  TECs  as  a  constant 
and  known  electric  current  is  applied  as  input  and  only  the 
temperature  gradients  need  to  be  determined. 

Several  models  for  TEG  modules  can  be  found  in  the  literature. 
Numerous  analytical  [6,9]  and  numerical  [10-13]  steady  state 
models  exist  for  simulation  of  TEGs.  However,  fully  coupled  and 
complete  transient  analyses  are  seldom  presented.  Transient  anal¬ 
ysis  of  TEG  subjected  to  a  load  change  from  steady  state  configura¬ 
tion  has  been  conducted  using  commercial  software  [14,15].  Mitrani 
et  al.  [16]  considered  temperature-dependent  material  properties 
and  the  effect  of  thermal  and  electrical  contacts.  Furthermore,  most 
models  exclude  the  Thomson  effects  [17,18]  and  their  influence  on 
apparent  Seebeck  coefficients.  Crane  [19]  showed  the  differences  in 
the  TEG  behavior  at  steady  state  and  during  transients  using  an 
uncoupled  thermoelectric  model  that  used  a  constant  current  value 
in  the  thermal  balance  equation  and  no  Thomson  effects. 
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Nomenclature 

t 

Time  coordinate 

X 

Spatial  coordinate 

L 

Length  of  p-n 

R 

Internal  resistance 

u 

Thickness  of  ceramic  hot  plate 

Ro 

Load  resistance 

l2 

Thickness  of  ceramic  cold  plate 

I 

Electric  current 

A 

Cross-sectional  area 

V 

Voltage 

As 

Convection  area  of  heat  sink 

a 

Heat  rate 

k 

Thermal  conductivity 

w 

Stored  energy 

P 

Material  density 

h 

Film  coefficient  of  heat  sink 

C 

Specific  heat  capacity 

a 

Seebeck  coefficient 

Subscripts 

P 

Thomson  coefficient 

P 

p-type  material 

7T 

Peltier  coefficient 

n 

n-type  material 

£ 

Electrical  resistivity 

1 

Material  of  hot  plate 

T 

Temperature 

2 

Material  of  cold  plate 

Too 

Ambient  temperature 

In  this  paper,  we  examined  the  behavior  of  a  TEG  when  exposed 
to  unknown  and  varying  heat  source  representative  of  energy 
harvesting  scenarios.  We  simulated  the  power  generation  charac¬ 
teristics  of  the  TEG  under  specified  temperature  profiles  that 
emulate  varying  heat  flux  applied  to  the  hot  side.  The  model 
parameters  for  a  commercial  TEG  element  were  characterized  by 
correlating  the  simulation  results  and  the  experimental  observa¬ 
tions.  We  show  that  the  developed  numerical  model  captures  the 
thermal  behavior  and  the  electrical  behavior  of  the  TEG  in  both 
open  and  closed-circuit  configurations. 


2.  Theoretical  modeling 

2.1.  Assumptions 


ApCppp 


—  k  A  d2^P  .  ^VrZ  t 

dF  -  p  p  9x2  +AJ  -V 


dTp 

dx 


(1) 


1  at 


0XZ  /in 


(2) 


where  Aj,  Cj,  pj ,  kj,  ej,  fij  and  7}  are  the  cross-sectional  area,  specific 
heat  capacity,  density,  thermal  conductivity,  electrical  resistivity, 
Thomson  coefficient  and  temperature  of  material  j,  respectively.  I  is 
the  electric  current  flowing  in  a  closed-circuit,  t  is  the  time  variable 
and  x  is  the  spatial  coordinate.  Subscripts  p  and  n  denote  the  two 
TEG  elements. 


The  basic  assumptions  made  include  neglecting  heat  losses  due 
to  radiation  and  transverse  convection.  The  flow  of  heat  and  current 
in  the  TEG  is  assumed  one-dimensional.  All  materials  are  assumed 
homogeneous.  Material  property  anisotropy  is  taken  into  account 
by  selecting  properties  relevant  along  the  dimension  selected  for 
the  analysis.  Thermal  conductivities,  electrical  resistivity,  and 
specific  heat  capacities  of  materials  are  assumed  constant  within 
the  operating  temperature  range.  While  the  numerical  methods 
used  did  not  require  any  of  the  above  assumptions,  these  simpli¬ 
fications  help  in  understanding  the  fundamental  mechanisms.  The 
reduction  in  the  dimensionality  of  parametric  space  enhances  the 
clarity  during  the  correlation  of  experimental  observations  and 
model  findings. 

2.2.  Governing  equations 

The  schematic  of  a  thermoelectric  generator  consisting  of  two 
thermoelectric  elements,  p-  and  n-type  is  shown  in  Fig.  1.  The 
simulation  dimension  is  divided  into  three  local  frames,  namely  {xi, 
x,  X2}.  The  x-coordinate  is  chosen  with  the  origin  at  the  bottom 
surface  of  the  TEG.  The  length  of  the  TEG  is  denoted  by  L.  x\  and  X2 
are  local  coordinates  in  the  hot  and  cold  plates  at  either  ends  of  the 
TEG  elements,  respectively.  The  power  at  each  junction  induced  by 
the  Seebeck  and  Peltier  effects  and  direction  of  thermal  conduction 
through  the  two  elements  is  shown  in  the  figure.  The  current,  J(t), 
flowing  through  a  load  resistor  indicates  that  the  device  is  oper¬ 
ating  in  power  generation  mode. 

Governing  equations  for  p-  and  n-thermoelectric  elements  are 
derived  from  the  energy  balance  equation  for  an  infinitesimal 
element,  dx,  as  given  in  Eqs.  (1)  and  (2)  [17,20]: 


hA. 


(r2'-r.)  f  T, 


m 


Cold  plate  T2(t,  x2) 


TL(t) 
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x2=L2 
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x2=o 
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Hot  plate  Tfi,  xj 

T.  Q{t)  f  T/(t) 
J(t)  Rr, 
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Fig.  1.  Schematic  of  thermoelectric  generators  with  hot  and  cold  plates. 
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In  contrast  to  an  earlier  model  [20],  we  considered  that  the 
thickness  of  the  hot  and  cold  plates  is  in  the  same  order  of 
magnitude  as  the  thermoelectric  elements.  The  temperature  field 
of  these  plates  is  also  solved  simultaneously  with  those  of  ther¬ 
moelectric  elements.  Governing  thermal  balance  equations  of  hot 
and  cold  plates  are  shown  in  Eq.  (3)  with  subscript  j  =  1  for  the  hot 
plate  and  2  for  the  cold  plate. 


PjCj 


07}  0^7} 

6t  J  dxf 


(3) 


2.3.  Initial  conditions 


The  initial  temperature  fields  of  elements  are  set  to  the  ambient 
temperature,  Too.  The  subscript  j  indicates  the  material  domain  with 
j  =  1,  2,  p  or  n  denoting  the  hot,  cold  and  the  two  TEG  segments, 
respectively. 

Tj(0,Xj)  =  Too  (4) 


2.4.  Boundary  and  interface  continuity  conditions 

The  boundary  condition  applied  at  the  bottom  surface,  X\  =  0, 
corresponds  to  the  applied  unsteady  heating  rate,  Q(t).  The  top 
surface,  *2  =  has  the  film  convection  boundary  condition 
employed.  At  the  interfaces  between  two  materials,  x\  =  L\  or  x  =  0 
and  x  =  L  or  X2  =  0,  both  the  temperature  and  heat  flux  are  assumed 
continuous  implying  a  perfect  bond  between  the  TEG  elements  and 


the  hot/cold  plates. 

X!  =  0  Q_(t)  =  (5) 

X!  =  Li(x  =  0)  T,(t,L,)  =  Tp(t,  0)  =  T0(t)  (6) 

x1=L1(x  =  0)  h(t,L,)  =  Tn(t,0)  (7) 

x1=L1(x  =  0)  -A,k^(t,L,) 

=  ccmT0(t)I(t)  -  (t,  0)  -  M„|  (t,  0)  (8) 

x2  =  0(x  =  L)  Tp(t,L)  =  T2(t,  0)  =  Tl(C)  (9) 


unit  cell,  respectively.  These  values  are  obtained  by  dividing  the 
plate  area  by  the  number  of  p-n  units  in  the  TEG. 


2.5.  Voltage  generation  in  TEG 

In  the  closed-circuit  mode,  the  output  voltage  generated  is 
computed  as  per  Eq.  (14): 

V=Vemf-IR  (14) 

with 


vernf  =  T0apn(T0)  -  TLapn(TL) 


(15) 


_  epL  snL 

R  ~  7^  + 


(16) 


The  open-circuit  voltage,  Vemf,  and  the  internal  resistance  of  the 
Seebeck  element  are  given  by  Eqs.  (15)  and  (16).  With  a  resistive 
load,  Ro,  electric  current  (/)  in  the  loop  can  be  determined  using  Eq. 
(17). 


j  _  ^emf 


(17) 


We  assume  that  Thomson  coefficients  of  the  p  and  n  type  of  TEG 
elements  are  both  independent  of  temperature  [21].  The  apparent 
Seebeck  coefficient  is  taken  to  be  related  to  the  Thomson  coeffi¬ 
cients  through  the  Kelvin  relationship  as  given  in  Eq.  (18). 


a](T)  =  af  +  n] :\n(T/T[ef)  (18) 

In  this  equation  ajef  is  the  Seebeck  coefficient  of  material  j  at 
reference  temperature  Tre f.  The  reference  temperature  considered 
is  the  room  temperature  (298  I<).  Substituting  Eq.  (18)  into  Eq.  (13), 
we  obtain  the  Seebeck  coefficient  of  device  at  temperature  T: 


«Pn (T)  =  (a]}ef  -  «pef)  +  (MP  -  (19) 

We  solve  the  transient  system  of  governing  equations,  Eqs. 
(1)— (3)  and  (15)-(17),  simultaneously  using  an  explicit  forward 
differential  approximation  in  time  and  a  central  difference  scheme 
for  evaluating  the  spatial  gradients.  Numerical  simulations  are 
performed  to  determine  the  temperature  and  electric  current  as 
a  function  of  time. 


x2  =  0(x  =  L)  Tn(t,L)  =  T2(t,0)  (10) 

x2  =  0(x  =  L)  apnTL(t)I(t)  +  kpAp^(t,L)  +  knAn^(t,L ,)) 

=  -A2k2g(t,0)  (11) 

x2  =  L2  -A2k2^(t,L2)  =  hAs(T2(t,L2 )  -  T„)  (12) 

with 

^pn  =  0Lp  —  0Ly\  (^3) 

In  Eqs.  (6)— (12),  T0  and  TL  are  the  temperatures  at  the  junctions 
x  =  0  and  x  =  L,  respectively,  h  and  As  are  the  convection  film 
coefficient  and  surface  area  of  the  top  plate,  respectively.  A\  and  A2 
are  the  cross-sectional  areas  of  hot  and  cold  plates  for  a  single  p-n 


2.6.  Convergence 

The  time  step  size  that  ensures  convergence  for  parabolic  PDEs 
is  as  follows[22]: 

max(At)  <  |^Ax2  (20) 

where  At  and  Ax  are  time  step  size  and  mesh  size,  respectively.  For 
the  TEG  elements,  this  time  step  resolves  to  approximately  28  ps  for 
the  mesh  size  chosen.  Fig.  2  shows  the  effect  of  spatial  mesh  dis¬ 
cretization  and  the  time  step  on  the  convergence  of  temperature 
fields.  Four  time  steps  and  mesh  densities  are  shown  in  the  figure.  A 
time  step  of  50  ps  and  spatial  grid  spacing  of  50  pm  are  selected  for 
further  analysis  is  sufficient  to  simulate  the  model  with  reliable 
accuracy  while  saving  computational  time.  For  this  grid  size,  the 
relative  error  is  less  than  0.1%  compared  to  the  finer  time  step  and 
mesh  sizes  (20  ps  and  30  pm). 
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Fig.  2.  Temperature  at  the  top  surface  of  the  cold  plate  computed  with  four  grid  sizes. 


3.  Experimental  analysis 

An  experimental  setup  was  designed  to  establish  parameters 
required  for  simulating  the  performance  of  the  TEG  and  to  validate 
the  model.  Fig.  3  shows  the  schematic  of  the  experiment  with  a  TEG 
(CUI  Inc.,  model  CP85438)  on  top  of  a  heat  source  (Corning  digital 
hot  plate,  model  PC400-D).  The  source  temperature  was  measured 
by  a  K-type  thermocouple  with  diameter  of  2  mm  placed  under¬ 
neath  the  TEG.  The  thermocouple  is  sealed  by  conductive  thermal 
grease  to  assure  uniform  heat  distribution  on  the  TEG  hot  surface. 
The  top  surface  of  the  TEG  is  cooled  by  natural  convection  with  a  no 
heat  sink.  The  TEG  has  ceramic  hot  and  cold  plates  made  of  alumina 
(AI2O3)  sandwiching  the  p-n  thermo-elements.  It  is  also  sealed  on 
all  the  sides  by  silicone  rubber.  The  temperature  of  the  top  surface 
(the  cold  side)  of  the  TEG,  T2,  is  detected  by  an  infrared  sensor 
mounted  above  the  heater.  Remote  sensing  with  the  infrared  sensor 
does  not  disturb  the  natural  convection  on  the  cold  surface. 

Initially,  the  system  is  at  room  temperature  (measured  as  298  K) 
and  the  hot  plate  was  switched  on  with  a  target  temperature  of 
388  K  (115  °C).  The  transient  temperatures  at  the  hot  and  cold 
surfaces  of  the  TEG  and  the  output  voltage  were  recorded  simul¬ 
taneously  into  a  digital  computer  using  a  data  acquisition  system 
(National  Instruments-6009).  The  data  acquisition  was  set  at  a  rate 
of  0.5  s  for  about  700  s  at  which  time  the  system  attains  a  steady 


state.  Fig.  3b  shows  the  schematic  of  various  elements  in  the 
thermal  conduction  chain.  The  p-  and  n-  thermo-elements, 
alumina  hot  and  cold  plates,  the  heater  base  plate  and  the 
contact  resistance  between  the  heater  base  plate  and  the  TEG  are  all 
modeled  for  proper  correlation  between  experimental  and  simu¬ 
lation  results. 

Experimental  data  are  shown  in  Fig.  4.  In  Fig.  4,  the  tempera¬ 
tures  measured  on  the  actual  thermoelectric  generator  and 
generated  power  (VI)  calculated  from  output  voltage  are  presented 
for  four  cases:  in  open-circuit  (represented  by  load  resistance,  Ro  = 
00 :  inf),  in  closed-circuit,  and  four  load  resistance  values:  R0  =  10, 
Ro  =  1.5,  Ro  =  0.5  and  R0  =  0.05  Q,  respectively.  Fig.  4a  shows  the 
temperature  at  the  hot  and  cold  sides  for  all  four  cases  (left  axis) 
and  the  temperature  difference  -  AT  (right  axis).  The  temperature 
difference  is  shown  only  for  one  case  (Ro  =  0.05  Q).  Due  to  the  low 
current  flowing  through  the  circuit  (the  maximum  value  is  about 
0.04  and  0.07  A  for  R0  =  1.5  and  Ro  =  0.05  Q,  respectively),  the 
Thomson  and  Joule  effects  are  not  identifiable  from  the  tempera¬ 
ture  fields.  The  electrical  output,  however,  is  strongly  coupled  to 
both  the  thermal  behavior  and  the  load  resistor  value.  The  gener¬ 
ated  power  increases  to  a  maximum  but  then  decreases  to  reach 
a  steady  state.  The  power  generation  is  highest  when  the  load 
resistance  value  is  close  to  the  internal  resistance,  which  is 
measured  to  be  about  1.5  Q.  This  observed  experimental  behavior  is 
in  agreement  with  trends  reported  in  the  earlier  literature. 

4.  Material  property  determination 

The  TEG  behavior  simulation  requires  several  material  proper¬ 
ties  to  be  estimated  for  thermal,  electrical  and  thermoelectrical 
behavior.  The  geometric  parameters  and  properties  of  ceramic  hot 
and  cold  plates  are  determined  from  either  the  manufacturer’s  data 
sheets  or  by  direct  measurements.  For  the  alumina  (AI2O3)  plates, 
the  properties  were  obtained  from  an  online  source  for  96% 
alumina.  The  properties  for  the  p-  and  n-  thermo-elements  were 
obtained  from  the  literature.  Geometric  parameters  and  material 
properties  used  in  the  model  are  given  in  Table  1  for  a  single  TEG 
cell. 

To  convert  data  obtained  from  a  TEG  with  several  cells  into  that 
for  a  single  cell,  the  number  of  p-n  pairs  counted  in  the  actual 
thermoelectric  generator  and  their  parallel/serial  network 
connections  are  determined.  The  voltage  results  from  the  model  are 
converted  to  the  complete  TEG  response  with  the  assumption  that 
all  p-n  cells  are  operating  under  the  same  conditions.  With  N  p-n 


Fig.  3.  (a)  Schematic  of  the  experimental  setup,  (b)  model  used. 
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Ms) 


Fig.  4.  (a)  Temperatures  at  both  sides  and  their  difference  and  (b)  output  power  ob¬ 
tained  from  the  TEG.  Input  heat  flux  is  shown  for  comparison. 


Table  1 

Properties  of  thermoelectric  elements  for  a  single  TEG  cell. 


Parameters 

Value 

Unit 

Source 

L 

2.2 

mm 

Measured 

U 

0.8 

mm 

Measured 

l2 

0.8 

mm 

Measured 

Ap 

2.25  x  10"6 

m2 

Measured 

An 

2.25  x  10"6 

m2 

Measured 

N 

127 

Measured 

As 

1.26  x  10"5 

m2 

Measured 

A i 

1.26  x  10-5 

m2 

Measured 

a2 

1.26  x  10-5 

m2 

Measured 

Pp 

7700 

kg/m3 

[20] 

Pn 

7700 

kg/m3 

[20] 

Pi 

3960 

kg/m3 

[23] 

P2 

3960 

kg/m3 

[23] 

h 

2.5 

J/m2  K 

[24] 

Pp 

45 

HV/K 

Estimated 

Pn 

-27 

HV/ K 

Estimated 

< 

25 

HV/ K 

Estimated 

< 

-32 

HV/K 

Estimated 

eP 

12 

pQm 

[26] 

£n 

10 

pQm 

[26] 

'<P 

1.8 

W/m  K 

[25] 

kn 

2.2 

W/m  K 

[25] 

*1 

18 

W/m  K 

[23] 

k2 

18 

W/m  K 

[23] 

cp 

200 

J/kg  K 

[20] 

Cn 

200 

J/kg  I< 

[20] 

Cl 

880 

J/kg  I< 

[23] 

C2 

880 

J/kg  K 

[23] 

cells  electrically  connected  in  series  in  the  TEG  (N  =  128  for  the  TEG 
considered  here),  the  relationship  for  electric  quantities  between 
the  actual  and  the  single  TEG  cell  are  in  Eqs.  (21)— (23): 


vs  -  "O'' 

(21) 

jact  _  jsingle 

(22) 

J?gct  =  NRs0ing,e 

(23) 

By  matching  the  transient  theoretical  results  with  the  experi¬ 
mental  data  in  open  electrical  circuit,  the  thermal  and  electrical 
properties  of  the  TEG  can  be  characterized.  After  the  characteriza¬ 
tion  process,  the  two  results  are  compared  in  closed-circuit  with 
different  load  resistance  and  in  different  applied  heat  flux  condi¬ 
tions.  In  the  first  subsection,  the  effect  of  Thomson  coefficients  is 
discussed. 

The  effect  of  Thomson  coefficients  on  the  temperature  and 
electric  power  output  predictions  is  shown  in  Fig.  5.  The  temper¬ 
ature  and  voltage  predictions  are  for  closed-loop  circuits  with  10  Q 
and  1.5  Cl  load  resistors.  Fig.  5a  and  b  shows  that  the  temperature 
predictions  are  not  significantly  affected  by  the  Thomson  effect 
term  but  the  output  voltage/current  is  considerably  impacted. 
Without  Thomson  effect,  the  predicted  voltage  is  lower  20%  than 
the  cases  with  the  presence  of  the  Thomson  coefficient.  Inclusion 
of  the  Thomson  effect  with  characterization  of  the  Thomson 


Fig.  5.  Comparison  results  with  and  without  Thomson  coefficients  for  (a)  R0  =  10  Q 
and  (b)  R0  =  1.5  G. 
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coefficient  is  critical  for  TEG  models  predicting  power  generation 
characteristics. 

4.1.  Thermal  behavior 

Since  the  heating  rate  of  the  hot  plate  used  in  the  analysis 
cannot  be  measured  with  considerable  accuracy  due  to  contact  and 
radiation  losses,  the  temperature  at  the  bottom  surface  of  the  TEG 
can  be  accurately  measured  by  a  thermocouple  and  translated  to 
the  model  as  a  time-dependent  input  condition.  The  thermocouple 
is  read  using  a  high  frequency  data  acquisition  system. 

The  thermal  characterization  is  conducted  under  open-circuit 
conditions  in  which  Vemf  is  measured  along  with  both  the  hot 
and  cold  temperatures.  In  this  case,  the  non-linear  coupling  terms 
are  negligible  in  the  equations  because  the  electric  current  in  the 
circuit  vanishes  and  the  system  consists  of  only  the  thermal  diffu- 
sivity,  Seebeck,  and  Thomson  terms.  Simulations  are  performed  by 
varying  the  conductivity  of  p-type  material  (kp)  from  0.9  to 
3.6  W/m  K  and  that  of  the  n-type  material  (fcn)  from  1.1  to 
4.4  W/m  K.  This  range  spanned  50-200%  of  the  reported  values  for 
the  conductivity  of  the  TEG  elements.  Fig.  6  shows  that  both  the 
transient  and  steady  state  behavior  of  the  TEG  being  tested  matches 
well  with  kp  =  1.8  W/m  I<  and  kn  =  2.2  W/m  K.  These  values  are 
typical  for  Bi2Te3  elements  and  such  values  are  often  reported  in  the 
literature.  Note  that  the  simulations  included  the  contact  resistance 
and  heating  plate  inertia.  Cheng  et  al.  [20]  identified  that  the 
specific  heat  capacities  of  elements  (Cp  and  Cn)  range  from  154  to 
210  J/kg  K  and  that  p-type  and  n-type  elements  have  similar 
properties.  In  this  analysis,  the  specific  heat  capacities  of  two 
materials  are  assumed  constant  and  equal  to  200  J/kg  K.  The  density 
value  for  the  TEG  elements  is  set  to  7700  kg/m3  based  on  the  data 
identified  from  the  literature. 

4.2.  Thermoelectric  coupling 

After  the  thermal  properties  are  determined,  the  open-circuit 
experiment  data  for  the  Vemf  Are  used  to  characterize  the  Seebeck 
and  Thomson  coefficients.  Fig.  7a  shows  different  voltage  curves  for 
different  values  of  Seebeck  coefficients  while  keeping  the  Thomson 
coefficients  constant.  The  absolutes  value  of  Seebeck  coefficients 
were  varied  in  the  range  from  20  to  100  pV/K  and  keeping  the 


Fig.  7.  Open-circuit  characterization  of  Seebeck  and  Thomson  coefficients. 


Thomson  coefficients  constant.  The  Seebeck  coefficients  are  esti¬ 
mated  based  on  correlation  of  open-circuit  voltage  predictions  and 
experimental  observations.  Fig.  7b  shows  the  comparison  of  the 
experimentally  observed  voltage  behavior  and  the  predicted 
behavior  when  Thomson  coefficients  are  varied.  The  Thomson 
coefficients  of  two  thermo-elements  are  varied  and  the  best 
voltage  characteristic  match  occurs  at  a£ef  =  2.5  x  10-5, 
=  -3.2  X  lO”5,  n p  =  4.5  X  1(T5  and  n„  =  -2.7  x  10“5  V/K. 
These  values  are  used  for  closed-circuit  behavior  prediction  and 
validation  of  the  simulation  results.  These  values  differ  from 
reference  values  in  reference  impurities,  effects  of  wire  leads  and 
contact  resistances.  Both  the  peak  voltage  and  the  steady  state 
values  of  the  theoretically  predicted  and  experimentally  observed 
data  are  correlated  with  the  set  of  material  property  parameters 
characterized  with  the  open-circuit  experiment. 


4.3.  Electrical  properties 

Next,  we  establish  that  the  electrical  resistivity  of  the  TEG 
element  matches  that  reported  for  Bi2Te3.  A  closed-circuit  test  was 
employed  in  order  to  monitor  the  resistivity  effect  in  output 
voltage.  We  pick  the  case  with  load  resistor  of  1.5  Q.  The  resistivity 
values  for  the  TEG  as  obtained  from  the  literature  are 
ep  =  1.2  x  10~5  Qm  and  en  =  1.0  x  10~5  Qm  [20].  Fig.  8  shows  the 
correlation  between  the  model  predictions  for  various  resistivity 
values  and  the  experimentally  generated  closed-circuit  voltage.  The 
voltage  response  is  seen  to  match  reasonably  with  the  resistivity 
values  obtained  from  the  literature.  These  values  are  adopted  for  all 


N.Q.  Nguyen,  K.V.  Pochiraju  /Applied  Thermal  Engineering  51  (2013)  1-9 


7 


further  analyses.  We  now  use  the  theoretical  model  under  applied 
heat  source  conditions  with  several  load  resistors  to  determine  if 
the  TEG  behavior  predicted  by  the  model  can  be  correlated  with  the 
experimental  behavior. 

5.  Correlation  of  simulations  with  experiments 

Fig.  9a-d  shows  both  the  experimentally  measured  data  and 
simulation  results  from  coupled  thermo-electric  model.  In  this 
correlation,  we  captured  the  time-dependent  heat  flux  applied  to 
the  TEG  into  the  simulation  by  imposing  the  temperature  profile 
seen  in  Fig.  9a  as  a  Dirichlet  boundary  condition  on  the  hot  side. 
The  match  of  temperatures  in  Fig.  9a  was  therefore  forced  and 
the  temperature  profiles  obtained  experimentally  and  through 
simulation  are  then  compared.  Fig.  9b  shows  very  good 


Time  (s) 


correlation  between  both  prediction  and  the  observations.  The 
temperature  difference  (AT)  is  also  seen  to  match  with  time  in 
Fig.  9c,  indicating  that  the  thermal  inertia  in  the  problem  has 
been  properly  modeled.  Fig.  9d  shows  a  comparison  between  the 
open-circuit  voltage  measured  and  that  predicted  by  the  simu¬ 
lation.  The  simulation  is  seen  to  capture  the  voltage  development 
trend  during  the  transient  as  well  as  the  steady  state  voltage 
output.  Quantitatively,  the  peak  voltage  predicted  differs  by 
about  20%. 

The  TEG  was  then  placed  in  closed-loop  for  power  generation 
mode  across  resistive  loads.  Three  load  resistors  were  tested.  Three 
load  resistors  of  10  Q,  0.5  Q  and  1.5  Q  were  selected.  The  current 
generated  by  a  TEG  is  expected  to  be  maximum  when  the  internal 
resistance  of  the  TEG  equals  the  external  load  resistance.  Fig.  10 
shows  the  simulated  and  experimental  temperature  and  voltage 
drop  characteristics  across  the  TEG  element  for  the  three  load 
resistance  cases.  The  results  show  that  predicted  and  the  experi¬ 
mental  transient  temperature  profiles  and  output  voltage  match 
closely  for  both  the  cases.  Slight  difference  in  output  voltage  with 
the  0.5  Q  resistance  is  observed  showing  an  effect  of  lead  wires  and 
contact  resistance  in  the  experiment.  Generally,  both  the  steady 
state  and  peak  voltage  values  correspond  very  well.  The  close 
match  of  the  peak  and  steady  state  values  indicate  accurate  char¬ 
acterization  of  both  the  thermal  conductivity  and  the  thermoelec¬ 
tric  properties  of  the  device.  Any  contact  resistance  values  are 
lumped  into  the  thermal  resistance  and  the  electrical  resistance  of 
the  materials.  The  voltage  generated  with  the  10  Q  resistor  is  much 
higher  than  that  of  the  0.5  Cl,  indicating  the  order  of  the  magnitude 
of  the  internal  resistance  of  the  system.  Therefore,  it  is  critical  to 
change  the  number  and  configuration  (number  of  units  in  parallel 
vs.  serial)  of  Seebeck  elements  to  match  the  electrical  resistance  of 
the  load  connected  to  the  TEG.  Another  experiment  was  carried  out 
with  the  same  TEG  but  at  a  higher  hot  plate  setting  temperature  of 
433  K  (160  °C)  and  a  load  resistance  of  1.5  Q.  Fig.  10c  shows  the  cold 
side  temperature  rise  with  time.  A  close  match  between  the 


Fig.  9.  Correlation  of  simulated  and  experimental  behaviors  of  the  TEG  element,  (a)  Hot  Side  temperature  profile  is  matched  between  the  simulation  and  experiment,  (b)  Cold  side 
temperature  profile  correlation,  (c)  Temperature  difference  between  hot  and  cold  sides,  (d)  Open-circuit  voltage  correlation. 
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Fig.  10.  Correlation  between  predictions  and  experimental  observations  (a)  R0  =  10  Q 
and  (b)  R0  =  0.5  fi  (c)  R0  =  1.5  Q. 


predictions  and  experiments  is  also  observed  for  both  temperature 
and  voltage  curves  in  all  cases. 


6.  Summary  and  concluding  remarks 

A  coupled  formulation  of  the  thermoelectric  behavior  of  a  TEG 
with  Seebeck,  Peltier,  Thomson  and  Joule  effects  included  in  both 
the  thermal  balance  and  the  voltage  predictions  is  developed  and 
used  to  determine  the  power  output  characteristics  of  a  TEG.  Heat 
input  represented  by  temperature  is  applied  to  one  junction  of  the 


TEG  while  the  other  junction  is  exposed  to  the  environment  under 
natural  convection  conditions.  These  conditions  are  typical  in 
an  energy  harvesting  application  in  which  waste  heat  may  be 
converted  into  a  useful  electrical  energy.  Finite-difference  tech¬ 
nique  is  used  to  compute  the  transient  temperature  and  voltage 
profiles  in  the  TEG  when  connected  to  a  load  resistor.  In  contrast  to 
a  thermoelectric  cooler,  for  this  problem,  both  the  current  flowing 
through  the  circuit  and  the  temperature  profiles  are  coupled  and 
time-dependent  and  hence  must  be  solved  simultaneously. 

Using  the  model,  we  determine  the  effect  of  the  Thomson 
coefficients  on  the  temperature  distributions  and  the  power 
output.  While  the  effect  is  insignificant  in  the  thermal  domain,  the 
power  generation  is  impacted  significantly  by  the  Thomson  term.  In 
the  present  study,  the  current  is  low  enough  that  the  impact  of  the 
Joule  heating  term  is  minimal. 

We  used  an  open-loop  experiment  in  which  the  TEG  was  sub¬ 
jected  to  a  known  temperature  profile;  several  parameters  for  the 
material  behavior  model  were  established.  The  parameters  whose 
ranges  established  from  literature  review  were  fine-tuned  with  the 
experimentation.  Model  predictions  for  closed-loop  power  gener¬ 
ation  scenarios  correlate  very  well  with  the  experimental  analyses. 

The  peak  voltage  generation  of  the  TEG  subjected  to  a  heat 
source  on  one  boundary  and  natural  convection  to  atmosphere  on 
the  other  is  not  at  the  steady  state.  The  steady  state  voltage  is  nearly 
half  of  the  peak  voltage  delivered  by  the  device.  Therefore,  the 
optimal  operating  point  for  these  devices  is  not  at  the  steady  state 
but  at  the  point  of  the  highest  temperature  differential.  In  a  power 
generation  scenario,  the  heat  input  to  the  device  must  be  controlled 
to  hold  the  temperature  differential  at  this  optimal  point.  Design  of 
such  a  controller  requires  accurately  modeling  the  system  behavior 
using  a  simulation  framework  as  outlined  in  this  paper. 

Acknowledgements 

Portions  of  this  work  were  conducted  with  support  from  US/ 
ARMY  ARDEC  Contract:  W15QKN-05-D0011-0051.  Authors  grate¬ 
fully  acknowledge  the  support. 


References 

[1]  S.B.  Riffat,  X.  Ma,  Thermoelectrics:  a  review  of  present  and  potential  appli¬ 
cations,  Applied  Thermal  Engineering  23  (8)  (2003)  913-935. 

[2]  S.S.L.  Chang,  Energy  Conversion,  Prentice-Hall,  Inc.,  Englewood  Cliffs,  NJ,  1963 
(Chapter  3). 

[3]  G.S.  Nolas,  J.  Sharp,  H.J.  Goldsmid,  Thermoelectrics:  Basic  Principles  and  New 
Materials  Developments  (2001). 

[4]  Y.  Yang,  X.J.  Wei,  J.  Liu,  Suitability  of  a  thermoelectric  power  generator  for 
implantable  medical  electronic  devices,  Journal  of  Physics  D:  Applied  Physics 
40  (18)  (2007)  5790-5800. 

[5]  C.  Wu,  Analysis  of  waste-heat  thermoelectric  power  generators,  Applied 
Thermal  Engineering  16  (1)  (1996)  63-69. 

[6]  Y.Y.  Hsiao,  W.C.  Chang,  S.L.  Chen,  A  mathematic  model  of  thermoelectric 
module  with  applications  on  waste  heat  recovery  from  automobile  engine, 
Energy  35  (3)  (2010)  1447-1454. 

[7]  X.  Gou,  H.  Xiao,  S.  Yang,  Modeling,  experimental  study  and  optimization  on 
low-temperature  waste  heat  thermoelectric  generator  system,  Applied  Energy 
87  (10)  (2010)  3131-3136. 

[8]  R.  Yang,  et  al.,  Transient  cooling  of  thermoelectric  coolers  and  its  applications 
for  microdevices,  Energy  Conversion  and  Management  46  (9,  10)  (2005) 
1407-1421. 

[9]  S.  Angrist,  Direct  Energy  Conversion,  fourth  ed.,  Allyn  and  Bacon,  Inc.,  Boston, 
MA,  1982  (Chapter  4). 

[10]  P.  Ziolkowski,  et  al.,  Estimation  of  thermoelectric  generator  performance  by 
finite  element  modeling,  Journal  of  Electronic  Materials  (2010)  1—10. 

[11]  L.  Chen,  et  al.,  Effect  of  heat  transfer  on  the  performance  of  thermoelectric 
generators,  International  Journal  of  Thermal  Sciences  41  (1)  (2002)  95-99. 

[12]  D.  Ebling,  et  al.,  Multiphysics  simulation  of  thermoelectric  systems  for 
comparison  with  experimental  device  performance,  Journal  of  Electronic 
Materials  38  (7)  (2009)  1456-1461. 

[13]  J.  Xie,  C.  Lee,  H.  Feng,  Design,  fabrication,  and  characterization  of  CMOS 
MEMS-based  thermoelectric  power  generators,  Journal  of  Micro¬ 
electromechanical  Systems  19  (2)  (2010)  317-324. 


N.Q,  Nguyen,  K.V.  Pochiraju  /Applied  Thermal  Engineering  51  (2013)  1-9 


9 


[14]  M.  Chen,  et  al.,  Transient  behavior  study  of  thermoelectric  generators  through 
an  electro-thermal  model  using  SPICE,  In:  International  Conference  on  Ther¬ 
moelectrics,  ICT,  Proceedings  (2006). 

[15]  E.E.  Antonova,  D.C.  Looman,  Finite  elements  for  thermoelectric  device  analysis 
in  ANSYS,  In:  International  Conference  on  Thermoelectrics,  ICT,  Proceedings 
(2005). 

[16]  D.  Mitrani,  et  al.,  One-dimensional  modeling  of  TE  devices  considering 
temperature-dependent  parameters  using  SPICE,  Microelectronics  Journal  40 
(9)  (2009)  1398-1405. 

[17]  P.E.  Gray,  The  Dynamic  Behavior  of  Thermoelectric  Devices,  John  Wiley  & 
Sons,  Inc.,  New  York,  1960  (Chapters  2-4). 

[18]  D.  Astrain,  et  al.,  Study  of  the  influence  of  heat  exchangers’  thermal  resis¬ 
tances  on  a  thermoelectric  generation  system,  Energy  35  (2)  (2010)  602—610. 

[19]  D.T.  Crane,  An  introduction  to  system  level  steady-state  and  transient 
modeling  and  optimization  of  high  power  density  thermoelectric  generator 
devices  made  of  segmented  thermoelectric  elements,  Journal  of  Electronic 
Materials  40  (5)  (2011)  561-569. 


[20]  C.H.  Cheng,  S.Y.  Huang,  T.C.  Cheng,  A  three-dimensional  theoretical  model 
for  predicting  transient  thermal  behavior  of  thermoelectric  coolers, 
International  Journal  of  Heat  and  Mass  Transfer  53  (9,  10)  (2010)  2001- 
2011. 

[21]  W.  Seifert,  M.  Ueltzen,  E.  Mueller,  One-dimensional  modelling  of  thermo¬ 
electric  cooling,  Physica  Status  Solidi  (A)  Applied  Research  194  (1)  (2002) 
277-290. 

[22]  Chapra,  S.C.  and  R.P.  Canale,  Numerical  Methods  for  Engineers,  fifth  ed..  2005, 
New  York:  McGraw-Hill. 

[23]  www.Matweb.com,  Materials  Properties  data. 

[24]  F.P.  Incropera,  D.P.  Dewitt,  Fundamentals  of  Heat  and  Mass  Transfer,  fourth 
ed.,  John  Wiley  &  Sons,  Inc.,  New  York,  1996. 

[25]  D.M.  Rowe,  CRC  Handbook  of  Thermoelectrics,  CRC:  Press,  USA,  1995, 
pp.  145-230. 

[26]  K.H.  Lee,  O.J.  Kim,  Analysis  on  the  cooling  performance  of  the  thermoelectric 
micro-cooler,  International  Journal  of  Heat  and  Mass  Transfer  50  (9,  10) 
(2007)  1982-1992. 


